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Abstract 

In these lectures, given at the NATO ASI at Windsor (2001), applications of the replicas non- 
linear sigma model to disordered systems are reviewed. A particular attention is given to two sets 
of issues. First, obtaining non-perturbative results in the replica limit is discussed, using as exam- 
ples (i) an oscillatory behaviour of the two-level correlation function and (ii) long-tail asymptotes 
of different mesoscopic distributions. Second, a new variant of the sigma model for interacting 
electrons in disordered normal and superconducting systems is presented, with demonstrating how 
to reduce it, under certain controlled approximations, to known "phase-only" actions, including 
that of the "dirty bosons" model. 

1 Introduction 

Starting from the seminal papers of Wegner [[jj and Efetov [0], a field-theoretical description based 
on the nonlinear a model (NLctM) has become one of the main analytical approaches to various 
problems in disordered electronic systems. The ensemble averaging over all configurations of 
disorder is performed either using bosonic ^ |[j or fermionic [p| re-replicated fields and taking 
the re — > limit in the results, or using supersymmetric fields [|| H], The main advantage of 
this approach lies in formulating the theory in terms of low lying excitations (diffusion modes) 
which greatly simplifies perturbative and renormalization group calculations, and also allows a 
non-perturbative treatment. 

The first application of this approach was a derivation [1-4] of the renormalization-group (RG) 
equations of the scaling theory |Q| of Anderson localization. For such a perturbative approach, 
generalized later for mesoscopic systems [5^], both the replica and the supersymmetric methods 
are equally justified: they ensure the cancellation of unphysical vacuum loops in a diagrammatic 
expansion. 

However, it soon became conventional wisdom that there exist two sets of problems, for each 
only one of these methods being applicable. The first set of problems can be schematically specified 
with the following TOE0 model: 

i i ij 

Here the second term represents schematically both lattice and disorder potential, while the last 
represents all possible two-particle interactions (including the BCS one where relevant). To include 
the interactions, at least at the perturbative and RG level, the original fermionic replica [|]] has 

'TOE stands for "Theory of Everything". In order to be really a TOE model in a condensed-matter context, the 
Hamiltonian ([l]) should also include spin terms, as later in this paper in relation to superconductivity. Naturally, 
either with or without spin terms, there is no hope for a rigorous approach to, let alone the exact solution of, this 
model for an arbitrary disorder/interaction strength. 
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been generalized by Finkelstein [Q], The interest in this approach has been greatly enhanced by 
the recent discovery |Q| of an apparent metal-insulator transition in 2D disordered systems in zero 
magnetic field. Naturally, the TOE model and its a model implementation covers a much wider 
variety of systems and phenomena than in the scope of the original NLcrM [1-4] describing only 
non-interacting disordered systems. 

After a few earlier attempts [9—11], it has been recently demonstrated [|l2| [□]] that the Keldysh 
technique jl4| ] provides a viable alternative to the replica approach for interacting systems. How- 
ever, the latter still remains one of the best available tools for consideration of interacting electrons 
in disordered systems. In the very least, it is clear that there is no simple way of applying the 
supersymmetry method to a many-particle fermionic system, so that the choice is between using 
Keldysh technique or replica trick. 

On the other hand, there exist another set of problems for disordered electron systems where 
meaningful results are essentially non-perturbative (or at least look like that - see section |]below). 
In the absence of interactions, many of them have been solved the help of the supersymmetry 
method (see, e.g., Refs. [15-20] out of many), while the viability of the replica approach was 
seriously questioned. 

The first, and arguably most famous, of such problems was solved by Efetov Jl5| | who used the 
supersymmetric NLcrM to derive the two-level correlation function (TLCF) in the universal ergodic 
regime for electrons in disordered metallic grains. The results proved to be identical to those for 
eigenvalue correlations [|TJ in random matrix theory (RMT), as had much earlier been conjectured 
by Gor'kov and Eliashberg [53]. This was the first microscopic derivation of essentially a non- 
perturbative result for noninteracting electrons in disordered media which has opened the way to 
numerous new results (again, see Refs. [15-20] out of many - this paper is not a proper venue to 
give a comprehensive list, let alone even a short review, of them). 

In many ways, this first non-trivial "supersymmetric" result jl5| ] seemed to be an excellent 
illustration of why the replica trick could only be used within a perturbative approach. For the 
(easiest) case of the unitary symmetry (one of Dyson's symmetry classes corresponding to the 
absence of the time-reversal symmetry, e.g. due to an external magnetic field) the irreducible TLCF 
is given by 

M») = -^> (2) 

where uj is the distance between two levels in units of S/w and 6 is the mean level spacing. This 
result is valid in the ergodic regime, i.e. for u> <C g where g ^S> 1 is the dimensionless conduc- 
tance. For u> 3> 1, the TLCF averaged over fast oscillations could be readily obtained from the 
standard diagrammatic techniques [ 23 1 with l/u> being the perturbation parameter (or from the 
perturbation theory in the framework of either the supersymmetric or replica NLcrM). However, 
the non-perturbative factor sin 2 w cannot be restored from the perturbation series. Since the replica 
trick is well justified only within the perturbative approach, it might seem rather hopeless to obtain 
the result (||) within the replica approach. And indeed, quite involved calculations by Verbaarschot 
and Zirnbauer p4[ have shown that a direct application of the replica trick (using either the bosonic 
or fermionic NLcrM) has apparently not reproduced the TLCF given by Eq. (||). 

In section |2| following Refs. [25-27], I will show that introducing a proper replica-symmetry 
breaking within the fermionic NLcrM leads to the non-perturbative result (g), albeit only in the 
asymptotic region u> ^> 1. This raises hope that one might eventually apply the replica approach 
for obtaining non-perturbative results for interacting electrons in disordered systems (see Ref. [ ^8| 
for the first steps in this direction). 

The replica-symmetry breaking is not the only way to get (at least seemingly) a non-perturbative 
result from the replica NLcrM. In section || I will briefly summarise how quite non-perturbative cal- 
culations [18-20] (that could only be performed within the supersymmetric method) have repro- 
duced part of the results for long-range asymptotics of different mesoscopic distribution obtained 
much earlier [29-32] by using a perturbative RG approach valid within any variant of the NLcrM. 

Then, I will describe in section |] how to include interactions into a modern variant of the replica 
NLcrM and will use this model in section || to describe different approaches to superconductor- 
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insulator transitions in thin disordered films. 



2 Replica-symmetry breaking: a way to non-perturbative re- 
sults? 

In this section we consider non-interacting electrons in disordered media - the problem described 
in a standard way by the first two terms in Hamiltonian ([l]) (this is equivalent to the Anderson 
model of disorder). The aim is to reproduce non-perturbative oscillatory contribution to the TLCF 
ofEq. (g) using the replica NLcrM. 

Before describing the derivation, it is worth specifying what is meant here by the replica- 
symmetry breaking. It is the standard fermionic NLcrM which does contain correct non-perturbative 
oscillatory contributions to the TLCF. The replica symmetry is not broken at the level of choos- 
ing a correct saddle-point approximation but rather fully exploited by involving a set of additional 
(in comparison to standard considerations [Q]) saddle-point submanifolds within the usual NLcrM. 
The symmetry is broken in a rather natural way when accounting for the contribution of these 
submanifolds to the TLCF. For the unitary case, such a procedure gives the exact results for all u>, 
although a real reason for this [ [33| ] may be incidental. There is no doubt, however, that it gives the 
asymptotically {oj ^S> 1) correct results for all the three Dyson's symmetry classes. 

Note an important analogy: in the supersymmetric NLcrM Andreev and Altshuler [|4j have 
reproduced the large-cj limit of Eq. (|j) by the saddle-point calculation of the supersymmetric inte- 
gral that involved an additional saddle-point which "breaks" the supersymmetry in a similar way. 
The same integral has been exactly calculated by Efetov [|l5|] without any breaking of the super- 
symmetry. 

This analogy will be extended by showing (with the help of a trick resembling the replica- 
symmetry breaking in the meaning explained above) that the n — replica limit of the exact 
integral representation of the TLCF obtained in Ref. [ [24] ] for the fermionic replica NLcrM leads to 
the same asymptotically correct oscillatory behavior of Eq. (Q). 

2.1 Formulation of the problem in terms of NLcrM 

Consider the two-level correlation function (TLCF) defined by 

R 2 {lu) = \{u{e + u)v(e)) - I. (3) 

Here (. . .) stands for the ensemble averaging, i.e. averaging over all the realizations of the impurity 
potential V^(r) in Hamiltonian (jl|), v(e) is the electronic density of states per unit volume defined 
in terms of the spectrum {e a } for a given sample as v(e) = L~ d J2 a <K £ — £ a), v = {v(e)) = 
1 / L d 5, and L, d are the sample size and dimensionality. All energies are measured in units of 6/n, 
in which the TLCF has the form of Eq. ( ]2j). In these units R 2 is expressed via the product of the 
retarded and advanced Green's functions as follows: 

iZaH^pfe&H-l], 5 2 (o;) = (G'"( e +|)G a ( £ -|)) . (4) 



The function S 2 can be expressed in the standard way [ |15[ p6[ | in terms of a generating functional 

S 2 {u;) = — lixq i~(Z B (u)), Z n =jv^e lS (5) 



n— >o n 2 du 

S[$M = /dr$(r)[|- (!+«?) A + V(r)]v>(r), V - +0 . (6) 

Here £ = p 2 /2m — Sp, ty a and ipa are conjugate 2n - component fermionic (Grassmannian) fields 
(corresponding to n replica components and 2 retarded-advanced components - further 2 spin 
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components are redundant in this section but will be included in section Q), A = diag(l, — l) ro ® 
11 n . [] The replica trick, n — ► 0, allows one to write the above expression for the product of Green's 
functions in terms of derivatives of Z (rather than of In Z) which makes the ensemble averaging 
straightforward. Performing it in the assumption that V(r) is a Gaussian <5-correlated random 
potential with the dispersion (27rz/T e i) _1 , where t c \ is the mean free time of elastic scattering from 
impurities, one obtains the following quartic term standing in the action S for ipVip: 



Sdi 



1 



dr [tp(r) • ip(r)] 



V>(r) 



(7) 



The next step is the standard Hubbard-Stratonovich decoupling that allows one to single out all 
"slow" modes of electron motion. In the absence of any symmetry breaking (due to an external 
magnetic field, or magnetic impurities, or spin-orbit coupling), there are two modes, equally con- 
tributing to quantum interference effects: the diffuson mode, corresponding to small transferred 
momenta, and the cooperon mode, corresponding to small sum momenta. It is instructive to per- 
form the decoupling directly in the r space, thus emphasizing that the diffuson mode corresponds 
to the normal pairing while the cooperon - to the anomalous pairing: 



exp{-S d i s } = 
Vc 



Vd 



exp 



dr 



4tv 



-tr e£(r) 



cxp < — / dr 



nl/ 4-1 \ , s 1 

^tr C +(r)c(r)- — 



2r i 

[VSc w (r)^- V>c+(r)V] 



(8) 



The diffuson field d u (r) is a matrix field of the same symmetry as ip ® ip while the cooperon field 
is of the same symmetry as ® ip. These two fields can be combined into the single 4N x AN 
matrix field cr(r) of the form 



c 

d T 



where d 



d, c 1 



(9) 



The extra 2x2 "diffuson-cooperon" space, explicit in Eq. (|9j), is usually referred to as the "charge 
conjugation" space (sometimes, as the time-reversal space): apart from being Hermitian, the matrix 
field a has the charge conjugation symmetry, 



= a = Ca*C 



* = (C*) J 



(10) 



where ^ = A= ('0 + , —ijj T ) , ^ T = -j= (tp T , ip + ), so that the second equation in ( flO| ) defines 
the charge-conjugation matrix C. The disorder functional (|8|) can be written in such an extended 
space as follows: 



exp{-S , dis } = 



exp < / dr 



8r 



-tr (7 



2r n 



(11) 



The functional integration measure in Eqs. (||) and ( pT[ ) is irrelevant in the replica limit. After 
performing the Gaussian functional integration over the fermionic fields in Eq. (^), one reduces the 
generating functional Z n to Z n = ~' F ^) with (• • -) CT standing for the functional integration 
over all the independent components of the field a, and the "free energy" is 



d d r 



8r o: 



-tr a 



— tr In 

2 



| + i& i A 



2r c i 



(12) 



A saddle-point condition for this functional, found by varying the field a is 



i~nvo = (r\ G \r) 



G 



(13) 



The same notation for A will be used in all extended spaces introduced later, assuming that it is a dyadic product 
of Pauli's a z in the RA space and of the unit matrices in all other sectors of the extended space. 
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In the class of spatially-homogeneous fields a, this condition is obviously satisfied by the field 
a = A as can be easily seen by transforming Eq. ( |l3| ) into the reciprocal momentum space and 
performing integration over all momenta in the r.h.s. of this equation in the pole approximation, 
justified provided that 

%T e l > 1 (14) 

What is usually not stressed is that the saddle-point condition (|l3]) is satisfied in the same approxi- 
mation not only by A but by any spatially-homogeneous diagonal matrix commuting with A with 
eigenvalues equal to ±1. For most applications (but not for what follows) the existence of this wide 
class of the saddle-point solutions at u> =^ is irrelevant. The reason is that for u = Eq. (fij ) is 
satisfied by a much wider class of matrices, a = Q, where Q satisfies the following constraints ^ 



Q 2 = ll 2 „, trQ = 0. (15) 

These constraints are resolved by representing Q as follows: 

Q = U^AU = 7^ AT, T = exp( °? t ' V (16) 



-it 0„ 

Since Q obeys the charge conjugation condition (|To[), U is a symplectic matrix (a unitary matrix 
whose elements are real quarternions [Q]), i.e. U <E Sp (2n). Then T is obtained by factorizing 
matrices U with respect to redundant matrices R E Sp (n) x Sp (n) that commute with A, i.e. U = 
RT, which reduces T to the form given in Eq. ( jl^ ) with t being an arbitrary nxn real-quarternionic 
matrix. This means p7| ] that Q belongs to the compact Grassmannian manifold (coset space), 
Sp (2n)/Sp (n) x Sp (n). This class of symmetry corresponds to Dyson's orthogonal class in 
RMT I will discuss other symmetry classes slightly later. 

The final step of the NLctM derivation is the expansion of tr In in Eq. ( p"2| ) in small u> and VQ 
around a spatially-homogeneous zero-frequency saddle point of Eqs. < JT~5| ) and (p"6]). To this end, 
one performs the similarity transformation in the tr ln-term representing it as 



tr In | Gq 1 - U[l C/ f ] + ^wlIAtf 



where Go is the Green's function of Eq. ( Jl 3[ > at uj = and a = A. After this, expanding the tr In 
to the first power in u> and the second power in the gradient operator £ = p 2 /2m — £p i=s t^n • V, 
performing the pole integration over ^, justified in the already used approximation £pT c i ^> 1, and 
neglecting in the same approximation tr a 2 term in Eq. (|l2"|), one finally arrives at the non-linear 
sigma model functional 



\ D{ VQY - ^AQ 



(17) 



where D — v F 2 T c \/d is the diffusion coefficient; at d = 2, the dimensionless conductance g — 
2tt 2 vD. 

If time reversal invariance is broken by a magnetic field or magnetic impurities, or spin-rotation 
symmetry is broken by the spin-orbit interaction, the above functional remains the same but the 
matrix U in Eq. ( |T6| ) belongs, respectively, to the unitary group U (n) (Dyson's unitary class) or 
to the orthogonal group, SO (n) (Dyson's symplectic class). Here a = 1 for the orthogonal class 
and a — 2 for the unitary and symplectic classes. This factor arises because unitary and symplectic 
classes have been obtained from orthogonal [Q| by the suppression of massive modes corresponding 
to the time-reversal or spin-rotational symmetry breaking and a subsequent reduction of the Q 
matrix rank. The coefficient a also absorbs an extra factor in the symplectic case due to the 
redefinition of the mean level spacing S in the chosen units. 



It is worth stressing that the condition tr Q = arises from choosing an equal number of replicas in both retarded 
and advanced Green's functions G r and G a . With this choice, non-zero values of trQ would correspond to 
"massive" modes that will not be contributing to the results. 
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For simplicity, I consider in the remaining part of this section only the unitary class, referring 
to the original publication J2(| ^] for the other two classes. I will also limit considerations to 
the ergodic regime corresponding to the level separations much smaller than the Thouless energy, 
E T ~ D/L 2 (which in the chosen units coincides, up to a numerical factor, with the dimensionless 
conductance g). In this regime the gradient term in Eq. ( |l7| ) may be neglected, and the NLcrM 
functional reduces to the zero-dimensional limit: [ [I~5| l 



-Tr [AQ] , 



(18) 



with Q becoming a spatially homogeneous matrix. 

2.2 Calculating the asymptotics of the TLCF 

Limiting all further considerations to the ergodic regime only, Z(uS) given by Eqs. (§) and (|T8j> can 
be represented as 



Z n (uj)= jVQcxp -z— Tr AQ 



(19) 



where the measure is defined by 



vq = Yl dft™ da 



ra * 
ij 



dn = dT ■ T~ 



(20) 



Here T is the matrix parameterizing Q, Eq. (|l6|), and r and a refer to the replica indices which 
originate from G r and G a , respectively. In the large-w limit this integral is mainly contributed by 
the extrema of the functional which obey the standard condition [A, Q) = 0. This condition is 
satisfied by any matrix of the form Q = diag(Q' , Q a ), where Q r and Q a are the n x n Hermitian 
matrices whose eigenvalues are ±1 and Tr (Q r + Q a ) = 0. This defines a highly degenerate 
saddle-point manifold which consists of G% n submanifolds specified by a particular distribution of 
n eigenvalues '+1' and n eigenvalues between Q r and Q a . These submanifolds can be divided 
into Ti + l classes of equivalence, Q p — diag(Qp, Q°), labeled by Tr Q r p = — Tr Q p = n — 2p, 
with p = 0, 1, . . . n. The p-th class has weight (CP) 2 , with CI = (™) . 

The matrix Q r p with (Qp 2 = 11„ and TrQ^ = n—2p can be parameterized by analogy with 
Eq. (@as 



T P = (Tp)^ XpTp , A p = diag(U Tl _ p ,-H p ), r; = ex P 



On— p 

-(r)t 



(21) 



where t r is an arbitrary p x (n—p) matrix. This defines the coset space G p = U (n) /U (n— p)xU (p), 
i.e. Therefore, Q p — diag(Qp, Q°) belongs to the manifold G p x G p and can be parameterized as 



Qp — TpA.pTp . 



T p -dia g (r;,T;) 



A, 



diag(A p , -A p ) 



(22) 



The integer p specifies the replica-symmetry breaking, as it describes the number of the —1 eigen- 
values in each Q r block (equal to the number of the +1 eigenvalues in each Q a block): in the 
symmetry-unbroken case, p = 0, and hence retarded and advanced blocks, Q r ' a , contain only 
positive or negative eigenvalues, respectively. 

Now one needs to take into account contributions from 'massive' modes (with mass oc l/u>, 
not to be confused with the massive modes oc neglected upon the derivation of the NLcrM) 

in the vicinity of each manifold (p2[). In the large-w limit these contributions may be considered as 
independent and the partition function is then represented by the sum of all of them: 



n 

Z B ( W ) = £(CS) S 

p=0 



T>Qexp —i—TrAQ p 



(23) 
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This expression is somewhat symbolic, as Q p = UA p U\ covers the entire symmetric manifold 
of the NLctM, Eq. (§J), including all the massive modes: indeed, TtUA p WA Tr A p T^AT = 
Tr A p Q, where we have substituted U = RT, as defined after Eq. (|l6|). This can be justified only 
as a perturbative (in 1/uj) procedure: a possible overlapping of massive modes originated from 
different manifolds is irrelevant in the large-o; limit, and each of the integrals in the sum (^3|) can 
be calculated independently of the others. 

Each term in the sum (^3j) contains both massive and massless modes. Indeed, we have used 
above the factorization U = RT with R being block-diagonal matrices commuting with A. The 
matrices T in Tr A p Q = Tr ATA p T^ still contain the subset of matrices commuting with A p that 
correspond to the massless modes. Therefore, we need to parameterize T in a way which enables 
us to factorize out these massless modes and perform the integration over the massive ones. 

The most suitable parameterization of T, analogous to that used in Ref. p-ij ] for the bosonic 
NLcrM, can be obtained by expanding the matrix exponent in Eq. (|l6|). By introducing matrix 
B = i(t t t)" 1 / 2 sin VtU, we represent T and thus Q = AT as follows: 

/ TUBBt) B \ _ ( H„ - 2BB^ BK(B*B) \ 

\ -flt K{B ] B) ) \ B^H{BB ] ) -(!„ - 2&B) ) ' K ) 



where 1Z(X) = \Jl n — X. The matrix B in this parameterization is not unconstrained, though. 
The Q — condition is fulfilled only when the matrices 1Z{BB^) and TZ(B'B) are Hermitian. 
This is so only when all the eigenvalues of BB^ and B^B do not exceed unity. Only under this 
constraint does Q, parameterized as in Eq.^4], still belong to the coset space U (2n)/U (n) x U (n). 
Nevertheless, this parameterization is very convenient. First, the corresponding Jacobian is equal 
to one so that the measure of integration ( |2Q| ) can be written as 

VQ = Y[ dB i:j dB*j = VB . (25) 

In addition, the representation of all the exponents in the sum (^3j) in terms of B is also very simple, 
Tr A p Q = 2(n - 2p) - 2Tr X p {BB^ + B^ B), so that we obtain: 

n 

Z n (u) = Y,( C n) 2 -^ (2p ~ n)z Z(")> (26) 
ZI{lo) = J VBexp [iluTt X p (BB^ + B^ B)] . (27) 

The region of integration in ( p7[) i s restricted by the constraint described after Eq. (|4|). Last, but 
not least, the parameterization^^) allows one to separate out the massless modes, which obey the 
condition [T, A p ] — 0, in each integral (p7|). Indeed, this condition is satisfied by all matrices T 
constructed from B which anticommute with X p , i.e. have the off-diagonal block structure. 

This means that in the representation of B in the block form reflecting the structure of X p = 
diag(l„- P) -Hp), 

the matrices B\,2 represent the massive modes, and 612 massless. When the massive modes are 
suppressed (Bi = and B2 = 0), the T matrices in Eq. ( p4[ ) constructed from p x (n—p) matrices 
61,2 only, parameterize the same degenerate p-the manifold, G p x G p described in Eq. (|22]), as one 
expects. 

By substituting the representation ( |28| ) into Eq. (p7[), we reduce to the product of integrals 
over the massive and massless modes: 



Zl{u) = j VB 1 VB 2 exp -2io;tr (^sj - B 2 b\) J Vb{Db 



(29) 
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Here the region of integration over b\^ depends on B12 due to the constraint on the eigenvalues of 
the matrices BB^ and B^B in the representation (^8|). Since the integral over B\2 is contributed 
only by the region where both tr B\B\ and tr B2B\ < l/o; <C 1, in the leading in 1/oj approx- 
imation we may put both B\ t 2 to in the constraint of the integration region over the massless 



modes b\^- In this approximation, as we have noticed after Eq. (28), matrices bi t 2 parameterize 
the p-th manifold ( j22| ) so that 

j vb = j Vb x Vb2 = J VQ P = r> 2 (G p ) , (30) 

where the measure of integration over T>Q p is defined in terms of T p in the same way that T>Q is 
defined in terms of T, Eq. (|2C|), and Q(G p ) is the volume of the compact coset space G p . This 
volume is expressed via the well-known volumes of the unitary group, il(U (n)), as follows: 



fi(U(ra-p))fi(U(p)) 

p 



7T \P 

(32) 



= (2 7 r)^" a -(^) 2 -^n r r r( lt i) -V (31) 

In the same large-w approximation, the variables B12 parameterizing the massive modes are un- 
constrained. Then the Gaussian integral over the 2[(n — p) 2 + p 2 ] independent massive modes 
yields 

// \( n— p) 2 
VBexp 2iwtr (BiB] - B 2 Bh =[— I 
L J icJ / Vicj/ 

Combining Eqs. ( f32| ) and (|T]) and omitting an irrelevant overall factor which goes to 1 when n — ► 0, 
we arrive at the following expression which is essentially the same as that derived in Ref. [ ^6| via 
the Itzykson-Zuber integral: 

W = E^ a - (2h>) (n-p W » ^ sC Sn r( n + 2-i) - (33) 

Here the summation over p has been extended to 00 since = for all integer n > p. This 
allows one to take the replica limit, n — * 0, in each of the terms in Eq. ([33|). Due to the fact that 
FP cx n p as n — > 0, only the terms with p = and p = 1 in Eq. ( p3| ) contribute to ^2 in Eq. (||). 

Let us stress that the replica symmetry is broken only now, in the n — > limit. Indeed, for any 
integer « ^ contributions of the terms with p and n — p are complex conjugate to each other, 
but for n — > we no longer treat them on equal footing. Thus, the result for Z n (oj) below is no 
longer a real function. Note, however, that in order to treat S2 (w) for all uj one should imply the 
u> — > u> + iS substitution which results in Z n (u>) being no longer a real function for any n. 

Omitting all the terms with p > 2, one obtains 

iujn piu(2 — n) 

Z >) = ^ + " 2 yMHl- (34) 

Substituting this into Eqs. (||) and (Q) and keeping the leading in 1/u terms only, one arrives at the 
expression ([|) for the TLCF. Although this expression coincides with the exact one (implying the 
above mentioned substitution to — > uj + iS), it has been actually derived only in the large-w limit, 
as is the case of the 'supersymmetry breaking' method of Andreev and Altshuler J34]]. 



2.3 The large-u; limit of the Verbaarschot-Zirnbauer Integral 

In order to obtain an explicit multiple-integral representation of the 'zero-mode' partition function, 
Eq. (|l9|), one uses the following 'polar' decomposition of Q: 

■+ ■ ' A vT^AV* 



Q 



VI - A 2 e~^ -A \ u 2 
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where — 1 < A^ < 1 and u\ 2 are unitary matrices. The appropriate measure of integration is given 
by 

DQ = Y[(\ z - \j) 2 JJdA; dcf>i dfi( Ul ) d/i(u 2 ) . 

i<j i 

The action in Eq. ( |l9| ) depends only on A, and the integrations over d/i(u) give the volumes of the 
appropriate unitary group. Thus one obtains 

+ 1 

Z„(w) = fi 2 (U (n)) I 6 2 (X) He-^dXt , S(X) = l[(X 3 - A,) . (35) 

This is equivalent (with accuracy up to factors going to 1 in the n — > limit) to the representation 
for S2 given in Eq. (2.24) of the paper by Verbaarschot and Zirnbauer, [ p4] ] which has been used 
for the critique of the replica trick. We will show that it leads, in the very least, to the exact large-cj 
assymptotic behavior of ^(w). 

The leading in l/u> contributions to this highly oscillatory integral (which does not have station- 
ary points inside integration region) come from the end points. To single out these contributions, 
we must take some A's close to +1 and the rest close to —1, which imitates replica symmetry 
breaking. Let us choose n—p of A's close to +1 and p of A's close to —1. Then we can split up the 
Vandermonde determinant in the following way: 

n 

5 2 {\) = Y[ \Xj - Ai| « 2 2p( - r ^A 2 + A 2 _ 

where 

n—p n 

A 2 + = l[ \X t - X,\ A 2 _= [] (36) 

— 1 i,j—n—p+l 

Reducing the integral (|5|) to the sum of such contributions, we represent it as 

n "hr n—p 

Z n {u) « fi 2 (U(n))^(C^)2 2p ("-P) / AlJJe-^dXj 

+°° 

X 

-1 



A 2 _ Y[ e'^dXj. (37) 

j—n—p+l 



Since in each of the integrals all the variables are close to one of the limits of integration, the 
second limit was extended to infinity. Now we make substitutions A, = 1 — Xi in the first integral, 
and A, = — 1 + x, in the second one, reducing the above sum to the form: 

n 

Z n {u) » il 2 (U (n)) J2(°n) 2 Mn - p) e^ 2 ^I n ^Ip (38) 
P =i 



where I p are integrals of Selberg's type: | |21 1 

J p 

Ip= A 2 J J d.r.v (39) 

i= x 

Substituting the known Selberg integrals and discarding an overall factor which goes to unity in 
the replica limit we arrive at 

w = Era 2 -^^, n=c*i[ ^ (40) 

p—0 j—1 v ' 
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This expression is exactly the same as Eq. 33 obtained in section 2.2 with the help of the replica- 
symmetry breaking. Therefore, the exact representation ( [35| ) does contain the true oscillatory 
asymptotic behavior of the TLCF. 

The authors of Ref. p4| ] have also drawn attention to the fact that there is an apparent contra- 
diction between the co = limit for 52 obtained from the replica trick and the exact supersymmetric 
result. Indeed, if uj is put to in the expression for S*2 following from Eqs. ( [35| ) and (||) the n — > 
limit is taken after that, one obtains 52 (oj = 0) = — 1. This cannot be correct as 5fte 52 (uj — > 0) > 
as follows from the definition (Q) and, moreover, it is known that S2(uj — > 0) — > <5(u). What is in- 
teresting, however, is that if one separates the singular, S^ 6 ^), and regular, S2 CS (uj), parts of the 
exact 82(0^), then S^iuJ — ► 0) = — 1. Therefore, the replica method gives 52 CS (cj = 0) correctly, 
and it is just 52 ing which is missing. However, the fact that Z n (uj — ► 0) is finite for any integer n 
does not necessarily implies that it is also finite (as a function of uj — ► 0) in the replica limit. For 
example, if the expansion of 52 (u>) before taking the n = limit contained a term proportional to 
uj , it would be singular in the replica limit. In other words, if a non-trivial dependence on the 
order of limits n — > and uj — > existed, in a spirit of the replica trick the n — > limit should 
be taken first. At the moment, though, this remains only a speculation. However, the fact that the 
large-w limit of this integral reproduces the correct results (Q) makes it plausible that there are only 
technical difficulties rather than one of principle in the application of the replica method. 

2.4 The small uj limit 

The purpose of this section was to demonstrate explicitly that non-perturbative oscillatory contri- 
butions to the TLCF of electrons in a random potential could be extracted from the standard NLerM 
formulated in fermionic replicas Q| many years ago. To this end, all one needs is to parameterize 
all the non-trivial saddle-point manifolds corresponding to the broken replica-symmetry and de- 
scribing 'massless modes' of the theory, and expand the action in the vicinity of these manifolds to 
include 'massive modes'. The very similar approach has been used in the supersymmetric NLerM: 
the non-perturbative oscillations have been extracted by the expansion around two extremal points 
one of which breaks the supersymmetry [p4]]. Since the exact supersymmetric calculation of the 
TLCF was well known (JT^j, it was clear that the supersymmetry breaking [|34j] was just a con- 
venient method of extracting the large-w limit (and going beyond the universal 'zero-mode' ap- 



proximation). It has also shown in section 2.3 that the exact integral representation of i?2(w) does 
contain the correct behavior in the large-w limit. 

The small uj limit is, arguably, more interesting as it governs a number of non-perturbative 
results obtained within the supersymmetric NLerM. I will briefly summarise some of them [18-20] 
in the following section, emphasizing that they have reproduced some of the results obtained much 
earlier [29-32] by using a perturbative RG approach valid within any variant of the NLerM. 



3 Tails of distribution functions 

All extensive physical characteristics of disordered systems fluctuate from sample to sample. How- 
ever, only in mid-eighties it was understood that the scale of such fluctuations is governed by the 
quantum coherence length so that at T — > they are not reduced with increasing the sample size. 
Immediately after the discovery p8| , [39| ] of the universal conductance fluctuations (UCF), it was 
shown [|9[ that their distribution function is Gaussian in its bulk, but has long lognormal tails. The 
existence of such tails appears to be a common feature of distribution functions of many observable 
quantities [29-32] such as global and local density of states, different relaxation times, etc. With 
increasing the disorder, the part of these tails in distribution functions was increasing so that the 
entire distribution of any local quantity was becoming lognormal in the critical regime in the vicin- 
ity of the metal-insulator transition [^lj]. The characteristic feature of the lognormal distribution 
is that the logarithm of its s th moment is proportional to s(s— 1) InL where L is the coherence 
length (at zero T, this is just a sample size). Therefore, all the moments (or in the limit of weak 
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disorder, where only the tails are lognormal, all the high moments) scale with different exponents, 
i.e. system shows multifractal behaviour. 

The results for the tails of the distribution functions that had been originally obtained[29-32] 
by the RG treatment of an extended (in a way described below) NLerM, have later been reproduced 
[18-20]in a much more elegant way in the framework of the standard supersymmetric a model. I 
do not intend to describe this approach here but only want to stress that some of the steps seem to 
be absolutely impossible within the replica treatment. The Q field in the supersymmetric model 
belongs to a "supermanifold" that includes both compact and a noncompact sectors. The existence 
of the noncompact sector is absolutely crucial for finding the distribution tails. With properly 
imposed boundary conditions required for finding the tails, Q = A does no longer represents the 
saddle-point. A correct saddle-point equation reduces, e.g., to the Liouville equation [|l9|] for a 
single scalar parameter that parameterises the (only relevant) noncompact sector of the space of 
the field Q (instead of, e.g., the parameterisation ( [24] ) in terms of the unconstrained matrix B or 
any other matrix parameterisation required in the replica a model). A non-perturbative solution to 
this equation provides for the existence of lognormal tails of the local density of states distribution 
[ fZ0| l or, equivalently, for the multifractality of the wave function [|l9|. Had such a solution been 
found before the RG solution of Refs. [29-32],it would be taken as yet another 'example' of the 
incapacity of the replica method. However, although the replica method works perfectly within 
the RG approach, there is still no clear understanding of why the results obtained within the non- 
perturbative approach outlined above are exactly the same as those within the RG approach to the 
extended NLerM (in the limit of the weak disorder, to which the validity of the non-perturbative 
supersymmetric approach is limited). 



3.1 Local density of states in open systems 

The two-level correlation function considered above is not well defined for u> < g for an open 
system, as level widths become of the order of Thouless energy, Et = TtyVerg which is much 
bigger than the level mean spacing 5 (their ratio is of order g 1 in the metallic limit). Its 
direct analog in this regime is the variance of the density of states (DoS) at the Fermi energy. It is 
useful to consider also the entire distribution of DoS which can be found in terms of its irreducible 
moments (cumulants). The relative values of the DoS cumulants can be expressed via the effective 
field-theoretical functional in a way similar to that in Eqs. (Q)-(^): 



_ (("*)) _ y (») 2S d 2S 

{v) s ~ n™ (2n 2 ) s do; 2 



The cumulants of local DoS are given by the same expression with partial derivatives substituted 
by variational ones, lu being considered in this case as a spatial-dependent source field in Eq. (||). 
As the higher derivatives with respect to li are involved in the above expression, the Tr In in 
Eq. (p~2() should now be expanded to the higher powers of to. This leads to the following additional 



contribution [B2fl to the er model functional ( 17 ): 



^ d [QM=- i i £r s Tr [.AQr +1 , ^=(^-) { j^§. (42) 

s— 1 ^ / \ /' 

As the bare value, T s q, of the additional "charge" is proportional to the s th power of the small pa- 
rameter T e i <5 ~ (Af)^ 1 ^/ L d , their direct contribution to the DoS cumulants is negligible compared 
to that obtained by the repeated differentiation (^T]) of the u> term in the standard a model.^| 



However, the charges T s sharply increase under RG transformations [32|. While "normal 



contribution to the DoS cumulants, i.e. those from the standard a model, Eq. (17), scale with the 



system size as appropriate powers of the DoS variance, the additional contribution from the vertices 



4 In calculating the local DoS cumulants, this small parameter is compensated by the appearance of powers of 
S(r — 0) — * nu/Td so that the contribution of the additional terms to the variance and higher cumulants of the 
local DoS is comparable [|3lh to that of the standard a model. 
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( |42| ) show the multifractal scaling due to this RG increase of T s . For Dyson's orthogonal class, it 
reads 

R s K .9o r «o e ^ ' , it = In — ss £ m — , (43) 

cr t 

where a is the conductivity at the scale of the system size L and cto is its bare value, that is the 
conductivity at the scale of the mean free path I. As the small parameter r s o scale as a linear power 
of s, for large enough s the e us ( s_1 ) factor becomes dominant. 

The additional contribution to the s th cumulant of the local DoS differs from that in Eq. ( |43| ) by 
the absence of the small parameter T s o. Such a contribution becomes dominant for s > u^ 1 (or for 
s > go in the weak disorder limit). The approximate equality in Eq. ( p3| ) refers to the weak disorder 



limit. In this limit, the result for the multifractal cumulants Q43J) has been reproduced within the 
supersymmetric approach [|l9j ^(J outlined above. Such a behaviour of the high cumulants leads 
to the lognormal tails of the distribution function. With increasing the disorder, u becomes of order 
1 in the vicinity of the metal-insulator transition in d > 2, or for the region of strong localization 
in d = 2 dimensions. In this region, the entire distribution becomes lognormal. The physical 
meaning of this is that there exist rare realisations of disorder which gives values of the local DoS 
(or, equivalently, local values of the wave function amplitude | W(r) | 2 ) which are much higher than 
the average value. The importance of this untypical realisations greatly increases with increasing 
the disorder so that they correspond to prelocalized states. 

3.2 Is there a paradox? 

The properties of tails of different distributions governed by the existence of the prelocalized states 
are described in detail in reviews [|2[ S3]. The point of this presentation is to underline certain 
peculiarities of the replica method. To this end, let me first outline how the RG results which 
lead to Eq. ( p3| ) have been obtained. The composite vertices (^) are not closed under the RG 
transformations of the a model given by Eqs. ( |l7| ) and (^2|). The index-permutation symmetry 
lacking in Eq. ( (42| ) is restored under RG transformations via generating the following additional 
vertices in each order: 



oo 



Tr (AQ)'] * , s t = s + 1 . (44) 

i=i l=i 

Thus in each s th order the RG equations become matrix equations, the matrix rank being equal 
to the number of partitions of the integer s into the sum of positive integers. These equations are 
exactly solvable [ plj ] for any integer n, and keeping only the largest eigenvalue leads in the replica 
limit n — ► to the result (|43|). As the RG approach is purely perturbative,| the applicability of the 
replica method is beyond doubt. However, there exist some apparent contradiction. 

The point is that one can choose different number of replicas for to represent the advanced 
and retarded sectors. In the simplest unitary case the Q matrix would belong to the coset space 
U (n)/U (n — to) x U (to) where n is the total number of replicas, and to is the difference in 
their numbers in the advanced and retarded sectors. In the replica limit, n — > and to — > 0, 
the result is naturally the same as before. However, the largest eigenvalue of the RG equations 
happens to depend only on n but not on to. Then, instead of taking the replica limit, one could 
fix to = 1 and take an arbitrary integer n. The resulting model is defined on the coset space 
SU (n)/U (n—1) isomorphic to sphere, i.e. in this limit the model becomes equivalent to the n- 



vector model. The composite vertices ( 42 ) reduce in this case simply to cos s ip where tp is the angle 



between the Ti-field and the direction of an external (magnetic) field. It is well known (see, e.g.. 



5 Eq. ( |43| ) has been obtained as a result of the one-loop RG analysis. The higher-loop contributions limit its validity 

— 3/2 

to s < g . Such a limitation ensures that no moment will be proportional to a positive power of the system 
volume, L d , i.e. all the multifractal dimensionalities remain positive. A similar restriction in the supersymmetric 
method of Ref. | fl^ ] follows directly from the inapplicability of the standard a model description at ballistic 
scales. 
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Ref. piQ) that such operators are irrelevant in the n-field theory, i.e. they only decrease under the 
RG transformations, in contrast to those of Eq. (|42|). So here is the contradiction as they should 
correspond, as described above, to a particular case of the matrix operators. It is clear, though, that 
a direct correspondence is anyway impossible, as the rank of the Q-matrices in Eq. ( p2| ) and thus 
the number of the eigenvalues in the RG equations depends on m and in the symmetric case m = 
it is equal to the number of partitions of the integer s into the sum of positive integers, i.e. much 
larger than in the case m= 1 when it is just equal to s. 

The resolution of the apparent paradox is in the fact that one should look only at correlation 
functions (observables) rather than at the RG dimensions of the composite vertices. In particular, 
only vertices [Tr AQ) S out of all the variety in Eq. (^4|) contribute to the cumulants (pj|), while 
all the rest serve to give correct eigenvalues. In calculating such a contribution, each eigenvalue 
enters with a coefficient of proportionality which is a polynomial of n and m. It appears that all the 
polynomials attached to "wrong" (in the case m = 1) eigenvalues vanish at all integer values of n 
thus restoring the correspondence between the matrix and the standard treatments of the n-model. 
Simultaneously, the n = m = replica limit gives the results (^3]) discussed above. On the face 
of it, it looks similar to the "replica symmetry breaking" discussed in the previous section: there 
exists a set of coefficients vanishing at all integer n and giving a nonzero contribution at n = 0. 
However, all the results discussed here have been obtained within the perturbative RG approach, 
and thus no trick was involved - the replica limit here is just to cancel unwanted vacuum loops in 
the RG diagrams. All these RG results can be equally easy obtained within the supersymmetric a 
model. 

One question that remains, though, is the following. The multifractal cumulants ( p3| ) have 



been obtained here within the extended a model, Eqs. ( |17[ ) and ( |42[ ), while later they have been 
reproduced within the standard SUSY a model. So, does one really need the additional vertices 

3.3 High gradient vertices: RG or not RG? 

The answer to this question is in getting - intermediately - even more vertices. Those in Eq. ( p^| ) 
have been obtained in expanding the Tr In term in Eq. ( |l2[ ) in higher powers of ui. But it could be 
also expanded in higher powers of VQ which yields the following set of vertices [|42|]: 

/2s 
Tr (j[d ai Q) d d r . (45) 
i=i 

The renormalization of these vertices results in the same increase of the largest eigenvalue as in 
Eq. (p3|), and actually, taken together with the w-vertex in the standard a model of Eq. (|l7|), they 
contribute to the cumulants, Eq. (p3|), in the same way (but for unknown - in any method - pre- 
exponential factors) as the composite vertices (^2]). Therefore, the latter are not even necessary to 
obtain the multifractality. 

An important point now is that, although the higher-gradient vertices (^5]) have been obtained 
from the higher-order expansion of the Tr In term in Eq. (|l2|), they are also generated in performing 
the RG transformation of the standard NLrrM. Usually they are omitted as being naively irrelevant. 
However, they become relevant at the one-loop level and actually govern long-tail asymptotics 
of different distributions. Note that in contrast to the composite vertices (ft4 ), the high-gradient 
operators remain relevant (in the same sense) also for the n vector model [[43||. 

Now the possible relation between the RG approach outlined here and the supersymmetric ap- 
proach of Refs. [|l9j ^(J is the following. The boundary conditions imposed to find correctly the 
tails of the distributions in the latter approach makes the saddle-point solution spatially inhomo- 
geneous, while in the former approach the relevance of the high gradients may indicate that such 
a spatial inhomogeneity develops spontaneously when starting from the standard homogeneous 
saddle-point. 

Finally, there is a possibility that the role of the high-gradient operators (^5]) is considerably 
more crucial. In the perturbative RG approach, they have only contributed to the higher moments 
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of conductance, DoS, etc, since the structure of the RG equations turns out to be triangular, i.e. 
the higher order gradient terms do not contribute to the renormalization of the lower order terms. 
However, in contrast to the composite operators (^2|), the operators ( |45| ) do not break any symme- 
try of the NLcrM. Therefore, one cannot exclude a possibility of a non-perturbative contribution 
from these operators to, say, the average conductance. Had such a contribution existed, it would 
totally ruin the one-parameter RG for the average conductance, and this the one-parameter scal- 
ing description of the Anderson localization. At the moment, such a possibility remains purely 
speculative. 



4 Coulomb and pairing interactions in the sigma model 

The main point of the previous sections is that the replica a model works perfectly for perturbative 
problems and - in certain cases - does reasonably well even for certain non-perturbative prob- 
lems within the non-interacting model. Since the inclusion of interactions beyond the mean-field 
approach is impossible within the supersymmetric method, the replica model the most reliable 
tool for models with interactions. The original fermionic replica a model [Q] has been general- 
ized to include the interactions by Finkelstein [Q]. It allowed him to reproduce earlier perturbative 
results (reviewed, e.g., in Ref. [[wj), to account where necessary for the Landau Fermi-liquid con- 
stants, and to derive the renormalization group equations describing (at least at a qualitative level) a 
metal-insulator transition in disordered interacting systems. This model has also been successfully 
extended [ p3[ | to allow for effects of Coulomb interaction in superconducting systems (like lower- 
ing the transition temperature T c by disorder - in an apparent, but well understood, deviation from 
the Anderson theorem [Q). The interest in this approach has been greatly enhanced by the recent 
discovery [||] of an apparent metal-insulator transition in 2D disordered systems in zero magnetic 
field. Although it is not at all clear whether the observed effects are, indeed, due the transition 
[ fi7[ and if they are - whether such a transition is, indeed, driven by interactions, a possibility 
of having the transition in a 2D disordered interacting system (schematically described by the TOE 
model (Q)), in contrast to a disordered noninteracting system, is intriguing by itself. Such a possi- 
bility is, undoubtedly, a driving force in a considerable revival of interest in Finkelstein's a model 
(see, e.g., Refs. jl^ , pjj |49[ |50|, |51 ]). At the moment, there is no clear evidence whether Keldysh 
techniques employed in Refs. [|l2f 13| | would give any edge over the replica method used through- 
out this chapter. Below I first describe how to include the Coulomb interaction (for simplicity, in 
a singlet channel only, as the inclusion of a triplet channel is technically almost identical) and the 
BSC pairing interaction into the derivation of NLcrM given in section |[ Then I will show (section 
Bh how to use such a model for describing the superconductor-insulator transition. 



4.1 Hubbard-Stratonovich decoupling 

Since the pair interaction does not conserve single-particle energy, in contrast to the elastic scatter- 
ing which was the only mechanism included up to now, the effective interaction functional should 
be dynamical. As usual, it is convenient to introduce imaginary time r implying everywhere the 
thermodynamic Gibbs averaging (together with the averaging over quenched disorder where appli- 
cable). Thus one considers thermodynamic Green's functions related (in the Matsubara frequency 
representation) to the retarded and advanced Green's functions of the previous sections by the stan- 
dard procedure of analytical continuation Then the effective functional corresponding to the 
interaction term in Eq. ([[]) can be written as 

<5>int = \ Jdxdx'ip s (x)ip s >(x')V xxl ^ sl (x')ip s (x) , (46) 

where x = r, r, and all the fermionic fields are anti-periodic in imaginary time r with period 1/T, 
s = (|, J.) is the spin index, V xx ' = S(t— T^)V c (r— r'), and V c represents the Coulomb interaction. 
After the replication, all the fields in Eq. ( |4q ) naturally have the same replica index. Similarly, the 
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BCS functional describing the attraction in the Cooper channel can be written as 



S sc = A I dx^{x)ipi{x)tl)i{x)^{x), 



(47) 



where Ao is the BCS coupling constant. 

The Hubbard-Stratonovich decoupling of the functional (^6|)-(^7|) is similar to that for the dis- 
order functional (0). To allow for all slow modes in the interaction functional(po|) one should in- 
troduce three matrix fields: the fields $ and / to take account of a small-angle scattering (a singlet 
channel) and large-angle scattering, and yet another one that corresponds to the Coulomb repul- 
sion in the Cooper channel. The last one would lead to the standard renormalization of the BCS 
attraction (see, e.g., Ref. [|5^]) and is not considered here; further, it is assumed that systems under 
considerations still have an effective attraction in the Cooper channel after such a renormalization. 
Under these conditions, the decoupling has the following form: 



-S in 



£>$exp<{-^ jdxdx' ®(x)V x J®(x') + i j dx ip a (x)®(x)ip„(x) 



-s sc _ 



T>f exp \ — — / dxdx' tr 



f(x)V x -}f(x') + i dx$ s (x)f ss ,(x)i> s ,(x) 



£>Aexpj-i J da;|A(a;)| 2 
i dx [A(x)i>t(x)i>±(x) - A(aj)^(a;)^ T (x)] 



(48) 



Now one doubles the number of components of the fields ip and ip as in Eq. ( 10 ) and performs the 
Gaussian integration to obtain 



S fi 



elds 



8r 



-Tra 2 
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-Tr In 



1 

i [ — a 

2r 



/ 



(49) 



Here the operator i equals if^ d T in imaginary time representation and becomes the diagonal matrix 
of fermionic Matsubara frequencies (s n = 7r(2n + 1)T) in frequency representation; Tr refers to 
a summation over all the matrix indices and to an integration over r and r (or summation over 
Matsubara frequencies in the frequency representation); S fields includes the ^-independent part of 
the action (48 ) quadratic in the fields $, / and A. The triplet channel, described by the field /, i s 
quite important: in particular, it can lead to the derealization in the presence of disorder f^, ^9|; 
however, such effects will not be considered here and this term will be ignored from now on as 
it is not relevant for the application of the model in section ^. The field a has the same structure 
as in Eqs. (§) and (|lo|), apart from explicitly including the 2x2 spin sector and replacing the 
2x2 retarded-advanced sector by the dependence on the imaginary time r (since it is diagonal 
in r, it becomes a matrix field in the Matsubara frequencies). In the absence of the fields $ and 
A, the saddle-point solution for the action (^) is formally the same as in the zero-temperature 
techniques of section || but the matrix A has a non-unit structure in the Matsubara (instead of 
advanced-retarded) sector: 



Q = C/ t At/, A = diag{sgne}. 



(50) 



The field $, corresponding to the singlet part of the Coulomb interaction, is diagonal in all the 
sectors. The "order-parameter" field A is Hermitian and self-charge-conjugate, diagonal in the 
replica indices and coordinates r and r, and has the following structure in the spin and time- 
reversal space: 



A{x) = \A(x)\e 



T2e 



FX(k)t3 



(51) 



where | A| and \ are the amplitude and the phase of the pairing field A(r, r), ti and f 4 sp are Pauli 
matrices (i = 0, 1, 2, 3 with f = 1) that span the charge-conjugate and spin sectors, respectively. 
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In the presence of the fields $ and A, the saddle-point equation for the effective functional ( \19\ ) 
can be formally written similarly to Eq. (Ob for the non-interacting zero-temperature case: 



im/<j(r) 



-£ + - — a + i[e + A + <S> 
2t c \ 



(52) 



Now a = A does no longer represent the saddle point solution for e ^ 0. Still, one can derive 
an effective functional by expanding the above Tr In within the manifold (^0|) in the symmetry- 
breaking field e+ A+$ and in gradients of Q, as has been done in the original works by Finkelstein 
[^J . An alternative is to make first a similarity transformation around A within the manifold ( pT)| ) 
to find the saddle point solution a sp in the presence of the fields. This can be formally done with 
the help of matrix Uq that diagonalizes the Hermitian field e + A + $: 



e + A + $ = Ul\U , a sp = W AU . 



(53) 



Here Uq belongs to the same symmetry group that defines the manifold (16) in the absence of the 
interaction fields. By substituting Eq. ( |53| ) into the saddle-point equation ( 52 ), one can easily verify 
that this is, indeed, a spatially-homogeneous solution, provided that Xt c \ ^> 1 which will always be 
the case for a dirty superconductor (Ar e i 3> 1) or a dirty metal (dimensionless conductance g ^> 1). 

Now one can perform the expansion of Tr In, Eq. (|49|), in gradients of Q and in the symmetry- 
breaking fields represented by the eigenvalues A and the matrix Uq, Eq. (|53|). It is convenient to 
employ the following parameterization: 



UlQU Q , Q = rfAU, 



(54) 



where Q represents the saddle-point manifold in the metallic phase and a is obtained from Q by the 
same rotation (|53|) as a sp is obtained from the metallic saddle point A. Therefore, Q is defined, as 
in the metallic phase, on the coset space S(2N) /S(N) ® S(N) where, depending on the symmetry, 
S represents the unitary, orthogonal or symplectic group as in the noninteracting case of section^. 
The parameterization (54) simplifies considerably all the subsequent derivations and leads to a 
new variant of the nonlinear a model which can be more convenient for many applications than 
the original Finkelstein's NLctM. After substituting a - 
obtains the following representation for the Tr In term: 



UqU^ AU U , Eq. (1541), into Eq. (|49|), one 



SS 



±Trln{C?o 1 + U [l W ] ~ i(U\U*)} 



where one can also include an external magnetic field with the vector potential A: 

- 1 - 1 

Z=7T- (p-eA) 2 -^. 
2m 



The expansion to the lowest powers of gradients and A is now straightforward and similar to that 
for the noninteracting case. It results after some calculations in the following effective action: 



S = i T r + J2 Jdr |A.| 2 + ^Tr 

The long derivative in Eq. (E^ is defined as 



| (dQ f - XQ 



A -ieAf 3 , Q 



d Q Q+[A ,Q] 



(55) 



(56) 



dQ = VQ 
where the matrix Ao is given by 

A = U d Ul , (57) 

and 8q = V — [ieAr3, . . .] is the long derivative ( |56| ) in the absence of the pairing field A. Both 
U Q and A should be found from the diagonalization of e + A + Eq. (53). Although such a 
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diagonalization cannot be done in general, it will be straightforward in many important limiting 
cases. For A = $ = 0, the field Ao vanishes, d —> do, and A — > e, so that the functional ( |55| ) goes 
over to that of the standard nonlinear a model for non-interacting electrons. 

The action ( |55| ) is most general in the current context. It allows one to develop a fully self- 
consistent approach to superconductivity of dirty metals in the presence of Coulomb interaction. 
However, any application requires a set of further simplifications. As a simple illustration, I will 
show below how to use the model in a dirty superconductor near the metal-superconductor transi- 
tion in the absence of Coulomb interaction. Then, I will show in section ^ how to use this model 
to describe the insulator-superconductor transition as a result of combined effects of disorder and 
interactions. 



4.2 Ginzburg-Landau Functional 

In the vicinity of the metal-superconductor transition one can expand the action ( |55| ) (in the absence 
of the Coulomb field $) in the pairing field A. A further simplification is possible in the weak 
disorder limit, g 1, when the Q-field can be integrated out to obtain an effective action for 
the A-field only. In the quadratic in A approximation, the kernel of this action will give, with 
due account for the disorder, an effective matrix propagator of the pairing field which governs 
properties of a disordered superconducting sample near the transition. 

To integrate over the Q-field, one splits the action (|5|) into S = So + 6>a where So is the stan- 
dard nonlinear a model functional as in the metallic phase. Then one makes a cumulant expansion, 
i.e. first expands e-( s o+ s &) in powers of Sa, then performs the functional averaging with e _5 ° 
(denoted below by (. . .)q) and finally re-exponentiates the results. The expansion involves only 
the first and second order cumulants since the higher order cumulants generate terms of higher 
order in A. Then the only contributions to the action quadratic in A are 



S cff [A] = ^£ |dr|AJ 2 -^(Tr(A- e )Q) ( 

-(^Tr[Ao,Q} 2 + ^(TrQd QA ^ ^ 

Expanding A and Ao to the lowest power in A and performing the functional averaging one finds 
[ pl] | the action quadratic in A as follows: 

5 eff [A] = |drA:(r)(r|/C w |r')A w (r'), (59) 



with the operator K,^ given by 



— -2nTY jlTr 

U e (u-e )<0 (. 



1 %e-u,\(°) C ' 



e(u)-e)< 



Heren^|(r,r') = (r|n c < rf |r') are the cooperon and diffuson propagators, with 



nf u , = [c + HJ , (6i) 

where C = —D (V — 2ieA) 2 defines the cooperon modes; fl d is obtained from fl c by putting the 
external vector potential A = 0. In the last term in Eq. (|50|), IljL (0) = nj^| (r, r); this term may be 
obtained by expanding in g~ x the cooperon propagator with the renormalized diffusion coefficient, 



C 



l- — nf„,(o) 
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which is a weak localisation correction to free cooperon propagator Ilf , (r, r') 
The summation over Matsubara frequencies in Eq. ( p0[ ) yields 



T 1 



V \2 AttT 



AttT 



(62) 



where To = T c q(B = 0) is the transition temperature of the clean superconductor in the absence of 
a magnetic field and ip is the digamma function. The weak localisation correction is proportional 
to the coefficient a u given by 



1 



1 



■kvV ^ Dq 2 
q 



Dq 2 



Dq 2 



4ttT 



AttT 



AttT 



For u> — the coefficient ao = a^^iT) can be simplified in the two limits: 



ao = < 



(63) 



f(l/2) 
8ir 2 vL d T ' 



where Tt = y/D/T is the thermal smearing length. 

The instability of the normal state (i.e. a transition into the superconducting state) occurs 
when the lowest eigenvalue of the operator K.^ becomes negative. The eigenfunctions of this 
operator coincide with the eigenfunctions of the cooperon operator C. The lowest eigenvalue of C is 
known to be Co = DB/cf>o, where </>o is the flux quanta. This ground state cooperon eigenfunction 
corresponds to the lowest eigenvalue JCo of the operator K. u . The condition /Co = implicitly 
defines the line T C (B) in the (T, _B)-plane where the transition occurs: 



In- 



To 
To 



Co 
AttT c 



aoCo 



(64) 



The term in the r.h.s. of Eq. ( p4| ) describes a l/g-correction to the main result. This weak lo- 
calisation is linear in the magnetic field B and vanishes as B — > as expected (the Anderson 
theorem [^6|). In a nonzero magnetic field the weak localisation correction to B c is positive which 
has a very simple explanation. The superconductivity is destroyed by the magnetic field when 
$(£) > $o, where is the flux over the area with linear size of the order of the coherence 
length £ ~ y/D/T and <l>o is the flux quanta. The weak localisation corrections reduce D and thus 
£. Therefore, one needs a stronger field to destroy the quantum coherence. The same reasoning 
explains the growth of T c in a fixed magnetic field. 

Note finally that it would be straightforward to include the leading weak-localization correc- 
tions in all orders of g^ 1 In by calculating the Q - averages in Eq. ( J58| ) via the renormalization 
group. This would lead to the renormalization of D in the cooperon propagator (^TJ), thus chang- 
ing the shape of the T C (B) curve. However, the value of T c (0) will again remain unaffected, since 
the superconducting instability is defined by the onset of the homogeneous zero mode in the oper- 
ator JC, Eq. (|62|), which does not depend on the value of the diffusion coefficient in the cooperon 
propagator. 

It is worth stressing that the Anderson theorem reflects certain properties of a model rather than 
those of real superconductors. If one allows for Coulomb interaction, then the critical temperature 
of the superconducting transition is no longer independent of disorder. Combined effects of the 
interaction and disorder lead to corrections to the transition temperature proportional (in a slightly 
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simplified way) to g" 1 ln 3 (T c r ol ) [g| ||, p. For a relatively weak disorder, the system still 
remains superconducting at T — > while for a sufficiently strong disorder the above corrections 
would suppress the superconducting pairing at any temperature and make a 2D system insulating 
at T — > 0. Such a mechanism of suppressing A by disorder gives a possible scenario for a widely 
observed superconductor-insulator (SI) transition in two-dimensional structures [55-58] which is 
most adequately described within the NLerM [[J5J ^| similar to those developed in this section. 
However, a widely accepted scenario for such a transition is based on very different, both in ori- 
gin and in implementation, models [59-63] where | A| remains finite and the superconductivity is 
suppressed by the loss of the phase coherence due to the fluctuations of the phase \ of the order 
parameter. I will show in the following section that these "phase-only" models can also be derived 
under certain parametrically controlled assumptions from the NLerM developed in this section. 



5 Superconductor-insulator transition in 2D systems: phase 
and amplitude fluctuations of the order parameter 

A wide variety of experimental data indicate the existence of a superconductor-insulator (SI) tran- 
sition in two-dimensional films [55-58] (see, in particular, [^8]] for earlier references). The tran- 
sition can be tuned by either disorder (changing with the thickness of a superconducting film) or 
magnetic field, thus being one of the most intensely studied examples of quantum phase transitions 

[0- 

One of the most accepted ways to understand the problem of the SI transition is based on the 
so-called Bose-Hubbard (or "dirty-boson") models [^3|. In this class of models exists a duality 
between charge-2e bosons (preformed Cooper pairs) and vortices. The superconducting phase is 
due to the bose-condensation of the charged bosons with localized vortices while the insulating 
phase is due to the bose-condensation of the vortices with the localized charged bosons. Another 
approach which captures the basic physics of granular superconductors is based on dissipative 
models |]65|] of resistively shunted charged Josephson arrays with the emphasis on the role of 
dissipation and Coulomb interaction. In both group of models [59-65], the transition is driven by 
fluctuations of the phase of the order parameter, i.e. the superconductivity is destroyed in spite of 
the existence of nonvanishing |A| locally. This approach seems to be rather different from that 
mentioned in the previous section where the disorder and interaction destroy the superconductivity 
by destroying | A| everywhere in a homogeneously disordered system. An experimental distinction 
between homogeneous and granular systems is not that strict[|55|] as it seemed a few years ago, and 
recent experimental observations [|5| ^] strongly suggest that the amplitude fluctuations in the 
vicinity of the SI transition are no less important than the phase fluctuations. 

The purpose of this section is to modify a general NLerM action for the description of granular 
systems with the BCS and Coulomb interactions. Such a description takes account of fluctuations 
of both amplitude and phase of the order parameter A, thus encompassing all the above described 
approaches. Both the Bose-Hubbard model [ |63| ] and the dissipative models [59-61] will be derived 
from this action via certain controlled simplifications made within the NLerM. The latter model is 
more general and allows one to go beyond certain limitations necessary in the derivation of the 
phase-only action. 

The starting point is a coarse-grained version of the TOE model Hamiltonian of Eq. (0) for a 
granular superconductor, where the kinetic energy includes terms iyajcij corresponding to tun- 
nelling hops between the grains. The derivation of the NLerM from such Hamiltonian follows the 
procedure described in the previous section. Thus a 'half-baked' effective functional obtained by 
integrating out the fermionic fields is similar to that in Eq. (^9|): 

S[a, A, $] = Tr cr 2 + — Tr | A| 2 + ^Tr ^f/" 1 * 
8r e i 4A 2 

(65) 

-i-i+-^-& + i(& + $ + e) . 



Tr In 

2 
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The difference here is that £ is the operator of the intra-grain kinetic energy (counted from the 
chemical potential), while t is the tunneling amplitude matrix (i.e. the inter-grain kinetic energy). 
Note that in the granular case the singlet field <£> is sufficient to decouple the Coulomb interaction 
term so that the triplet field / does not enter Eq. (|5|). All the bosonic fields have the same structure 
as before but for the addition of the m X m grain sector. The symbol Tr refers both to a summation 
over all these matrix indices and to an integration over intra-grain position r and the imaginary 
time t. 

The principal simplification for granular systems is that all the fields are spatially homogeneous 
inside each grain when the grains are zero-dimensional, i.e. their sizes L < £. Lt (£ and Lt 
are the superconducting and thermal coherence lengths) which is equivalent to |A|,T < l/r er g. 
For a diffusive grain r org = L 2 /D = H/Et while for a ballistic (chaotic) grain r crg ~ L/v ¥ . 
The tunneling matrix t = {Uj} depends only on grain indices, and the Coulomb interaction thus 
reduces to the capacitance matrix, U^ 1 — > dj/e 2 . 

To follow the procedure described in the previous section one needs, first, to solve (formally) 
the saddle point conditions, Eqs. (j53|) and (|4|) separately for each grain. As all the fields are 
spatially homogeneous inside each grain, Uo = Uo(i) commutes with the operator £. Then one 
only needs to expand the Trln to the first nonvanishing orders in and Aj, this expansion being 
justified when \t\, |A|,T <C 1/t c \ <C £ f . Thus one arrives at the following effective action: 

' • a 



S[0 ,A,*]=/ dT {£JM_ + £gi ti% } 



(66) 



where Sij = U (z)C/q (j), all the fields depend on t, Tr refers to all indices except i, j numerating 
grains, Si is mean level spacing in the i-th grain, and the tunneling conductance is defined by 
gjj = 2ir 2 \tij\ 2 /SiSj (whichis nonzero only for neighboring grains). Both A,; and Ua(i) (and thus 
Sij) should be found from the diagonalization procedure in Eq. (|53|), while Q is given by Eq. (p4|). 
The next step is to represent Uo(i) as 

U (i) = ^ e -3»(T)f 3 _ (67) 

This is similar to the gauge transformation suggested in Refs. p^ ] and used in Ref. [^7| to gauge 
out the Coulomb field. However, one cannot gauge out two inde pen dent fields, A and <£>. Substi- 
tuting the transformation ( |67| ) into the diagonalization condition (p3|), one reduces it to 

e + ^+A° = V+X t V t , (68) 

where A° is the field ( |5l| ) taken at \ = and the field $ is given in the r representation by 

= $ 4 - \d T Xi ■ 

Both $ and | A| are massive fields whose fluctuations are strongly suppressed. It is straightfor- 
ward to show that the fluctuations of $ are of order 8 which is much smaller than both T and |A|. 
Therefore, in the mean field approximation in $ this field can be neglected, $ = 0. This condition 
is nothing more than the Josephson relation in imaginary time [|60||. This locks the fluctuations of 
the Coulomb field <f> with the phase-fluctuations of the pairing field A, 

* = \d rX , (69) 

thus reducing the action (^6|) to one depending only on the fields Q and A. 

The mean field approximation in |A^| is valid for |A;| ^S> 5. It is equivalent to the standard 
self-consistency equation which formally follows from the variation of the action (pq) with respect 
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to |A|. In this approximation one finds |Ai| to be independent of i and r. Thus the first term in Eq. 
(pq) becomes a trivial constant, so that the action depends on Q and \ only: 



(70) 



The field \ m this action obeys the boundary condition x(t + /3) = x( T ) mod 2ir. In calculating 
the partition function, one should in general allow for different topological sectors corresponding 
to different winding numbers in x- 

The "phase-only" action (|70| ) includes neither fluctuations of | A|, nor fluctuations of $ beyond 
the Josephson relation, Eq. (p9|). Below this action will be reduced to the AES action. Still, we 
stress that the action ( J70| ) is more general than the AES action. Thus, in the absence of supercon- 
ductivity, A = 0, it was shown [ |&7| ] that the former contains a correct screening of the Coulomb 
interaction at low T, in contrast to the latter. This may also be important in the case when A is 
much smaller than the charging energy. 

To further simplify the action ( [70[ ) we note that the diagonalization conditions (|68|), in the 
absence of the | A | fluctuations, are the same for each grain and reduced to those solved in Ref . : 



V eE > 



)f 2 sin sgne 8 E - e 



(71) 



A = diag \/ e 2 + |A| 2 sgn e , cos ( 



V^ + IAI 2 



Then SV, in Eq. (tfOh can be expressed in terms of V as 



Si . 



Xij =Xi~ Xj ■ 



(72) 



Finally note that large- |e| contributions to the action ( |13| ) are strongly suppressed, while for |e| <C 
| A | one has A = | A|A which suppresses fluctuations of Q in each grain imposing Q — A. Then, 
all matrices in the action ([70]) are diagonal in the replica indices so that these indices become 
redundant. It finally reduces the action ( fToj ) to that depending only on one scalar bosonic field 
Xi(r), the phase of the order parameter, which is a function of only the imaginary time and the 
grain number. Indeed, the second term in Eq. (|70b reduces to a trivial constant; evaluating the 



(73) 



tunneling term with the help of Eqs. (17 lh and (|72h, one obtains 



where Xij (r) = Xi(-r)-Xj{r), X i} - 2 



2ffy JdrJdT'gliT-T 1 ) cos%y + gl (r-r') cos X 
Xij{r)±Xij{r') 



o o 

Uj — 2 

functions g n a (integrated over all momenta) are given by 

e sin et 



and the normal and anomalous Green's 



fJ 



\A 2 + |A| 2 



| A | cos er 



(74) 



The action (^3|) is exactly the AES action introduced in Refs. [ |59[ |60[ |6l| ]. Further simplifications 
are possible in two limiting cases. 

First, in the normal case ( A = one has in Eq. ( f74j ) 



g a = o. 



9l(r) 



rp2 

sin 2 ttTt 



(75) 



21 



Then the field x should be substituted, according to Eq. (j59|), by 2f r dr'$(r'). This corresponds 
to using the action ([73|), ([75]) for a set of normal tunnel junctions [p8[], as has been recently shown 
in Ref. [|67||; the functional ([70]) in the limit A = is equivalent to that of Ref. [|67|], Including 
disorder-induced fluctuations (i.e. going beyond the Q = A approximation) allows one to obtain 
[ [57t a correct low-T limit for the phase correlation function missing in the action (|73|). 

The action ( f70| ) is more general than that considered in Ref. $5%: although under the mode 
locking condition (^) it depends only on the fields x an d Q> me matrix Sy, Eqs. ( f7l| ) and (]72|), 
reduces to a simple f7 (1) gauge transformation as in in Ref. [ |&7| | only in the limit A = 0. 

The second limiting case, T <C | A | , is just the limit relevant in the context of the SI transition in 
granular superconductors. For T = 0, the summation in Eq. ( |74| ) can be substituted by integration 
which yields 

9n (r) = M Jfi(|A|r), g a (r) = ^ K (\A\r) 
ir 7r 

This is also a good approximation for a low-T case; substituting this into Eq. ([73]) gives the action 
for the dissipative model [|6C|, |s"H ]. Note that for |r — r'| <C |A| _1 , the main contribution in the 
tunneling action ( f73| ) is given by the normal term with the corresponding kernel proportional to 
|t — t'|~ 2 . The Fourier transform of this would give a term of the Caldeira-Leggett type[p5|[ 
proportional to \lu\. 

The tunneling action ( J73[ ) is non-local in r. As has been noted in Ref. [|o|| for the case of one 
Josephson junction, for sufficiently large capacitance the phase Xij changes slowly in comparison 
with | A| _1 , and in the adiabatic approximation x( T ') is changed by x(t) + (t'— t)9 t x(t). Making 
such an expansion, one obtains from Eq. ([73]) the following local action: 



S[x\ = / dr IE W*Ki ~ I cos XH } , (76) 

o « 

where Xi = d r Xi an d 



l C',; g% 3 + cos Xij 



u- 1 



4e 2 ^ |A| 8 

j 

Cij u) , 3 + cos Xij 



13 " 4e 2 |A| 8 

If all the self-capacitances are equal to C with E c oc e 2 /C being the charging energy, and all 
g\j — g*", then uu = J7 has the meaning of the renormalized charging energy. Ignoring a weak 
dependence of u on cos Xij m tne above relations, one obtains the renormalized charging energy: 

E c 

U = . (77) 

Here the coefficient # depends on the number of next neighbors for each grain, etc. A similar 
renormalization takes place for the next-neighbor off-diagonal energy . Now one can see that 
on the face of it the adiabatic approximation employed to obtain Eq. ([76]) is valid for U <§C |A|. 
However, in the region 3> | A|/£" c , where the charging energy ( f77| ) is strongly renormalized, the 
instanton-like solutions [69] may be important. This may further reduce the region of applicability 
for the local in r action (76). 

Finally, by introducing the operator n canonically conjugate to the phase x> one finds the 
Hamiltonian that corresponds to the action (f76|): 

& = E ~ \ A \9ij cos (Xi ~Xj) ■ (78) 
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This is just the Hamiltonian of the Bose-Hubbard model Q63Q which was first microscopically 
derived by Efetov [ f70| ] in the context of granular superconductors. 

To conclude, in this section the effective er model-type action ( |66| ) has been derived for a gran- 
ular system with zero-dimensional grains in the presence of the Coulomb and pairing interactions. 
This is the most general (in the present context) action which takes into account fluctuations of both 
amplitude and phase of the order parameter A. Neglecting fluctuations of | A| and fluctuations of 
the Coulomb field $ beyond the Josephson relation ( 69 ) reduces this action to the "phase-only" 
action (|70|) which still contains intra-granular disorder important for the correct screening for | A| 
small compared to the charging energy. Neglecting this disorder further reduces the action ( f7(i| ) to 
that of the AES model, Eq. (f73[). When the renormalized charging energy, Eq. ([77]), is much smaller 
than |A|, the action (|73|) finally goes over to that of the Bose-Hubbard model, Eq. (j78j), which is 
widely used for the description of the superconductor-insulator transition [|53|l . However, the above 
estimations show that this reduction is parametrically justified only for the region E c <C | A| where 
the transition happens at g t <C -E C /|A| <C 1 which corresponds to a strongly granular system. 
Note finally that the most general (in the present context) action ( |66[ ) describes both amplitude and 
phase fluctuations of the order parameter, being still considerably different from the NLcrM action 
for homogeneous systems. 
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